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This paper studies the stability of warm inflationary solutions when the viscous pressure is taken into 
account. The latter is a very natural and physically motivated ingredient of warm inflation and is seen to 
widen the stability range of warm inflation. The spectral index parameters, n^, ut, and their ratio are 
^^ , derived. The corresponding WMAP7 data are used to fix some parameters of the model. Two specific 
\^ • examples are discussed in detail: (i) a potential given by V{(j),T) = vi{<j)) + V2{T), and (ii) a potential of 
1^ ' the form ¥{(/), T) = avi{(f>) V2{T). In both cases, the viscosity has little impact on the said ratio. 
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I. INTRODUCTION 

Nowadays it is widely admitted that our Universe experienced a very early stage of accelerated 
expansion (inflation) which served, among other things, to produce the seeds that, in the course 
of the subsequent eras of radiation and matter dominance, developed into the cosmic structures 

iking, we might say there 
J], that ends in a very 



(galaxies and clusters thereof) that we observe today [l|, I4]. Broadly spea 
are two main competing scenarios of slow roll inflation: cold inflation [3 



i, 



that does not necessitate the latter 



short and violent phase of reheating, and warm inflation 
phase as the inflaton field decays into radiation at a sufficient (damping) rate, say F, during the 
accelerated expansion. This essential feature enables the Universe to smoothly proceed into the 
radiation phase, indispensable for the big bang nucleosynthesis. A further characteristic of warm 
inflation is that the radiation temperature satisfies T > H. As a consequence, thermal fluctuations 
result bigger than quantum fluctuations whereby the former likely are the main origin of cosmic 
structures [7[. 

Any inflationary model -whether "cold" or "warm"- must fulfill the requirement of stability; 
that is to say, its inflationary solutions ought to be attractors in the solution space of the relevant 
cosmological solutions. It means, in practice, that the scalar field, (p, must approach an asymptotic 
attractor characterized by ~ —{dV/dcj)) (2>H)~^ in cold inflation, and (j) ~ —(dV/d(j)){3H + T)^^ 
in warm inflation (see e.g. |8J, |9|). This ensures that the system will stay sufficiently near to the 
slow-roll solution for many Hubble times. Here V denotes the scalar field potential and H the 
Hubble expansion rate. 

In the case of warm inflation the conditions for stability have been considered by de Oliveira 
and Ramos [lO| and, recently, more fully by Moss and Xiong [11] who allowed the scalar potential 
and the damping rate to depend not only on the inflaton field but on the temperature of the 
radiation gas as well. This automatically introduces two further slow-roll parameters and renders 
the conditions for a successful warm inflationary scenario even less restrictive. 

The aim of this paper is to take a further step in the analysis of the stability of warm inflationary 
solutions by considering the presence of massive particles and fields in the radiation fluid as well 
as the existence of a the viscous pressure, H, associated to the resulting mixture of heavy and 
light particles. This novel quantity is well motivated on physical grounds. It can represent either 
a genuine viscous pressure (as in radiative fluids [l2;ll3|]), or an effective negative pressure linked 
to the production of particles by the inflaton field, or the decay of the heavy fields into light fields 



(see, e.g. 



14 



171]). or both. Whatever the case, it may significantly infiuence the dynamics of 



warm inflation. In particular, the inflationary region occupies a wider section of the phase space 



associated to the autonomous system of differential equations when 11 7^ than otherwise \ii 
and its effect on the matter power spectrum (of a few per cent) will be possibly detected in future 
experiments [19i]. Moreover, the adiabatic index 7 of the mixture will not necessary be 4/3 but 
some value in the range 1 < 7 < 2, that will depend on the proportion between heavy and light 
particles in the mixture. 

This paper is organized as follows. Section II presents the basic equations for warm inflation 
with viscosity. Section III analyzes the stability for inflationary solutions. Section IV determines 
the power spectrum produced by the thermal fluctuations. Corrections to the spectral index are 
considered when the viscosity comes into play. Section V applies the full theory to two speciflc 
examples differing in the type of potential but with identical forms for the viscous pressure and 
damping rate, T. The seven year WMAP data are used to constrain the model. Finally, Section 
VI summarizes and discusses our flndings. Hence forward a coma means partial derivative. We use 
units in which c = h = 8ttG = 1. 

II. BASIC EQUATIONS 

We consider a spatially flat, homogeneous and isotropic universe dominated by an scalar fleld (p 
and radiation, of temperature T, endowed with a viscous pressure 11 = —3 C,H, where C is the semi- 
positive deflnite coefficient of bulk viscosity. A potential V{(f),T) and the damping rate T{(f),T) 
dictate the evolution of the field so that the corresponding Klein-Gordon equation assumes the 
form 

4>+i3H + T),j) + V^ = 0. (1) 

In this scenario, the Friedmann equation reads 

3H^ =p^ + p^, (2) 

where p^ = (l/2)(/)^ + V{(I),T) and p^ = Ts, denote the energy densities of the scalar fleld (the 
inflaton) and radiation, respectively. Here 

s ~ -Vt , (3) 

is the entropy density. This expression follows from the thermodynamical relation s = —df/dT , 
when the Helmholtz free-energy, / = (1/2)0^ + V{(j),T) + p^ — Ts, is dominated by the potential. 
Since for an expanding Universe the viscous pressure is necessarily negative the total pressure 

p = i</)2 - Vi^, T) + (7 - l)Ts + H, (4) 



is lower than in the absence of viscosity. Obviously, the immediate effect of 11 is to further accelerate 
the expansion. Clearly, the adiabatic coefficient, 7, must slowly vary with expansion but, for 
simplicity sake, we will consider it constant. 

The conservation of the stress-energy can be expressed as 

Ts + 3H{-fTs + n) = r(/)2 . (5) 

Making u = (j) and introducing the dimensionless strength dissipation parameter, Q = T/{3H), the 
slow roll equations take the form 

"= 3ff(l + Q) ' ^^ = J ' 3F=y(0,r). (6) 

The number of e- folds is given by 

where (/>* and (/"end denote the value of the inflaton field when the perturbations cross the horizon 
(i.e., when k = aH) and at the end of inflation, respectively. 
Finally, in warm inflation the customary slow roll parameters, 

I fVA^ V.ss . V,sT., 






are supplemented by two others, 






that gauge the temperature dependence of the potential and the damping rate, respectively. 

III. STABILITY ANALYSIS 

We are interested in finding the impact of the viscous pressure, 11, on the stability of the slow 
roll. This immediately implies that we must focus on the strong regime (i.e., Q ^ 1) as in the 
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weak regime {Q <^ 1) viscosity will be practically non-existent 

To find the conditions for the validity of the slow roll approximation, we perform a linear stability 
analysis to see whether the system remains close to the slow roll solution for many Hubble times. 
In cold inflationary scenario, the slow roll equation is of first order in the time derivative. Choosing 
the infiaton field as independent variable, the conservation equations ([1]) and ([5]) can be written as 
first order equations in the derivative with respect to <j), indicated by a prime, 

x' = F{x) , (10) 



where 

/ u 

(11) 

V ^ 

Thus, the system ([1]), ([5]) becomes 

u' = -3H - r - Vc^u-^ , (12) 

s' = -SHjsu-^ - 3HU{Tuy^ + T-'Tu . (13) 

Here the Hubble rate and temperature are determined by ([2]) and ([3]), respectively. 

Taking a background x which satisfies the slow roll equations ^, the linearized perturbations 
satisfy 

dx' = M{x)5x - x' , (14) 

where 

(ab\ 

M=\ ], (15) 

is the matrix of first derivatives of F evaluated at the slow roll solution. Linear stability demands 
that its determinant be positive and its trace negative. 
The matrix elements read, 



Hi Q_ 

s I '^ [l + Qf 



B = -{-cQ- rHh^e + h{l + Q)\ , (17) 



Hs (^ e \ r _J]_ / 6(l + Q)^-2e 

[1 + QY)\ +7^7 V6(l + Q)2-e 



c=iii[^-Tr^T^]\^+^ r:..^j... : )r m 



^ H f ^ Qe \ HU ( Qe SH 1 

uV 7^(1 + Q)^/ -"7^7 1 7^1 + <3) 27^p^J 

In the strong regime (Q S> 1), the determinant and trace of M assume the comparatively simple 
expressions 

detM = ^2^f4-26 + c+(c-26)^-|^V (20) 

u^ V 7^7 2 7^4/ 



and 



trM = -S-3Q + 7(c - 4) + -^ 



9 n 

^7 



Sufficient conditions for stability are that M varies slowly and that 

6> 0, 



, , 4-3aV2 

c < — ^ 26, 



l + a 



(21) 



(22) 



where a = -3 — . Upon these conditions the determinant results positive and the trace negative, 

implying stability of the corresponding solution. Expression ()22[ 1) generalizes Eq. (27) of Moss 

and Xiong 11 1. 

Since the chosen background is not an exact solution of the complete set of equations, the forcing 

term in equation (J14p depends on x' , and will be valid only if x' is small. The size of x' depends 

on the quantities u/{Hu) and s/{Hs). From the time derivative of (IHS) we obtain 

H _ e 

lP~~l + Q' 

Combining this with the other slow-roll equations, dUl) and ((612), we get 



(23) 



u 

mi 



1 

A 



c[A{l + Q) - BQ] - 4 



and 



s _ 3 



1 + Q 

A{3 + Q) - B{1 + Q) 
1 + Q 



e + 



13 + {Ac - A)ri ^^ 4^6 



1 + Q 



e + 



Q 



1 + Q 



Aj3 - 2Ar] 



1 -/ 



{1 + Q)[Ac{Q - I) + Q + l]c 
(i-/)Q 



(24) 



where 



A = 4(l + Q) + ylc(Q-l), A 

3(1 + Q)2 C 



P7 + 7~'n 
p-y — KH 



B 



n 



Py — KH 



f 



c, 



K = Py- 



P-t 



(25) 

(26) 

(27) 



2 Q -/He' ' ' C 

Notice that when 11 — )• 0, one has that ^ — )• 1, i? — )■ 0, / — )• 0, and therefore the equations 
()23p - (l25p reduce to the corresponding expressions in [11]. Obviously, the value of the parameter k 
in this limit depends on the specific expression of the viscosity coefficient, (^; but it does not alter 
the value of B in the said limit. In this limit, the k parameter could take any value depending of 
the model. Its value does not affect the 11 —t- limit. 



IV. DENSITY FLUCTUATIONS 



The thermal fluctuations produce a power spectrum of scalar density fluctuations of the form 



m 



V, 



2 n2 



^/lTQ. 



(28) 



Note that the power spectrum of fluctuations in inflationary models where the friction coefficient 
depends also on the temperature, i.e., T = T{(j),T), was considered recently in Ref.[20(]. 
We calculate the spectral index by means of 

ns-l = —^ . (29) 

HVs 



By virtue of the equations (|23|) - (j25|) . we obtain 



Pie + p2(3 + psT] + Pib 
ns-l = , (30) 



where the pi coefficients are given by 



10(2 + Q) - A{3 + 5c + Q) + B {1 + Q + (5c/2)Q) 

^^ = ttq ' ^''^ 

AiQ - 1) - log ^3^^ 

P3 = ^'^ ' ^' ^^': ; - — ^ ' --'^ (33) 



(34) 





P2- - 


1 + Q 


8(1 


+ Q) 


- A{2 + 2C+2Q + 3cQ) 






1 + Q 


P4 = 


3(1 + 


Q)[l + (l + 5c/2)Q] 

/ -I n\ j^^ 



For Q ^ 1, and assuming c of order unity, the pi coefficients reduce to 



pi = -10 + A-B{l + 5c/2), p2 = A-W, p^ = 8-A{2 + 3c), p^ = ^^^^±^^ , (35) 

and A = (5(4 + Ac). Therefore (|30l) becomes 

_ _ 10 - A + i?(l + 5c/2) _ 10-^ 8-A{2 + 3c) 3(1 + 5c/2) ^ 

""' (4 + ^c)Q ' (4 + Ac)g^+ (4 + Ac)g "^^ (4 + 4c)(l-/)''- ^""^^ 



Ui, i.e., 



The tensor modes happen to be the same as in the cold inflationary models 

Vt = H\ (37) 

and the corresponding spectral index is 

„,-l = -^.. (38) 

With the help of of (j37p . (j28p and ([6jl) the tensor-to-scalar amplitude ratio can be written as 

, = ^^^(^'^) (39) 

9V^i/3T(l + Q)5/2- ^ ^ 

The recent WMAP seven-year results imply the upper-bound r < 0.36 (95% CL) [21] on the scalar- 
tensor ratio. Below, shall make use of this bound to set constraints on the parameters of our 
models. 



V. EXAMPLES 

In this section we apply the formahsm of above to two specific cases sharing the same damping 
and bulk viscosity coefficient but differing in the potential. In general, the damping coefficient may 
be written as 



.T) 



rp\ n 

^0. 



(40) 



with n and m real numbers and 0o, "^Oi ^ind Fq some nonnegative constants. The damping term 



23t |. From 



has a generic form given approximately by F ~ g^(/)^r, where g is the coupling constant 
Ref. 2j] the damping term, r = r((/), T), is related to the relaxation time of the radiation and for the 
models with an intermediate particle decay, r = t{(J)) is linked to the lifetime of the intermediate 
particle. Different choices of n and m have been adopted. For instance the case n = m = was 
considered by Taylor and Berera Ml, whereas the choice m = 2, n = —1 corresponds to the damping 



term first calculated by Hosoya 241]. This expression slightly differs from those in 
single index rather than two was considered. 

As for the bulk viscosity coefficient we use the general expression 



23 



25l |. where a 



(41) 



where Co is a positive semi-definite constant and A an integer that rnay take any of the two values: 



1/2 



n 



A = 1/2, i.e., ( oc p.y [2g] (see also Ref. J27l| ) and A = 1, i.e., C oc P7 uM 



A. First example: a separable potential 



We begin by considering the comparatively easy separable potential [11] 



VicP,T)=viiT) +t;2(</') 



where vi{T) = -^g,T^ - ^m^T^ and V2i^) = ImU^. 



(42) 



In this case, the strength of dissipation parameter and the slow rolls parameters are 

r(</>,T) 



Q 



y/3(f 1 + V2) 



(43) 



and 



4a2 






/3 



mm J 



V 



m. 



V 



6 = 0, c = n 



(44) 



respectively. 



Recalling that p^ = —TVt, we obtain the following expression for the radiation energy density 



mA 2 mO, 



(45) 



To get a lower limit to the friction term in the presence of viscosity, we follow the method of 
We cast Eq. ()28p in the form 

rni tt3 

and use ([6l2) together with (l40]l to arrive to the cubic equation 



m- 



(46) 



4X -52^-53 = 0. 



(47) 



with 



92 



12CP. 



9i 



45 Q* 



where Q* = Q yjl + Q and X = H/T. 

The only real solution of the Eq. ()47p is found to be 

X = Z sinh {9) , 

with 



Now, the condition for warm inflation, T > H (i.e., X < 1), translates into 

3 / 45 \ C " 



Q > g^Vs 



1 



(48) 



(49) 



(50) 



(51) 



2 \7r'^jg* J T^ _ 

Note that, as C ^ 0, the viscous pressure augments the range of Q, and that the condition for warm 
inflation without viscosity, Q > g^:Vs, is trivially recovered in the limit C — ^ 0. 
Using (HH) and (|15|) the inequality ^UT\i can be written as 



Q>g*Vs - . .. 

7 V 45 

Since Q > an upper bound on C,q follows, namely. 



^_ 3Co / 27rV y ^^4A-3 



l-A 



SU ^ Smax — q \ Ar 



(52) 



(53) 
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FIG. 1: Top row of panels: Plot of the tensor-scalar ratio r as a function of the spectral index ris, for two 
values of the A parameter in the case of example 1 (i.e., potential (02])). Bottom row: Same as the top row 
but assuming no viscosity (Co = 0). In each panel the 68% and 95% confidence levels set by seven-year 
WMAP experiment are shown. The latter places severe limits on the tensor-scalar ratio [21]. 



Figure [U depicts the dependence of the tensor-scalar ratio, r, on the spectral index, n^, for the 
model given by Eqs. ()40p . (j4ip . and ()42p when A = 0.5 and when A = 1. From Ref. [21], two- 
dimensional marginalized constraints (68% and 95% confidence levels) on inflationary parameters 
r and n.o, the spectral index of fluctuations, defined at /cq = 0.002 Mpc"-^. The seven-year WMAP 



data 2l| places stronger bounds on r than the five-year WMAP data 23] ■ In order to write down 
values that relate n^ and r, we used Eqs. (i36]l and (f39|) . and the values g* = 100, 7 = 1.5, 



Co = (2/3)ci^L, and m^ = 0.75 x 10"^, T = 2.5 x 10"^, Fq 



for rn, = 0,n = 0; and m^ = 2.5 x 10 ^, T 
(j)Q = 0.6 for m = 2, n = 0, in Planck units 



1.75 X 10~ 



M- 



1.2 X 10-^ TO = 3.73 X 10"^ (t)Q = 0.3 
Fo = 3.58 X 10-^ To = 5.63 x 10-^ 



Figur^D suggests that the pair of indices {m = 2, n = 0), corresponding to the right panel, is 
preferred over the other pair of indices {m = n = 0), left panel. Likewise, it shows that there is 
little difference between choosing A = 1 or A = 0.5 as well as with the case of no viscosity, i.e., 
Co = 0. 
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Panel in Fig. [U 


N 


r 


N 


r 


top left (to = n = 0) 


60 


0.351 


75 


0.314 


top right (m — 2, n — 0) 


60 


0.094 


75 


0.074 



TABLE I: Results from first example with A = 1 (The results for A = 1/2 are very similar). Rows from top 
to bottom refers to panels of Fig. [T]from left to right. 



Panel in Fig. [T] 


N 


r 


N 


r 


bottom left (to ~ n ^ 0) 


60 


0.350 


75 


0.318 


bottom right (to = 2, n = 0) 


60 


0.094 


75 


0.074 



TABLE II: Results from first example with no viscosity, i.e., Co — 0. 

Table [11] indicates the value of the ratio r for A = 1 and different choices of the pair of indices m 
and n when the number of e-folds is 60 and when it is 75. Very similar values (not shown) follow 
for A = 0.5. All of them can be checked with the help of Eqs. (j36p and (j39p . 

A comparison of the results shown in both Tables indicates that only in the case of the pair 
(ttt, = n = 0) with N = 75 (top and bottom left panels in Fig. [1]) viscosity makes a non-negligible 
impact. 



B. Second example: effective Coleman- Weinberg potential 



As a second example we consider the^ffective Coleman Weinberg potential near the critical 
point, for T < Tc in the region about (J)b 



where vi = k — T^log((5/T) with 5 ~ m-p\, and V2 = ( 
condition k > -^ = kmin must be imposed. 
The strength parameter adopts the form 

r(c/>,T) 



V{ct>,T) = avi{T)v2{(t>), (54) 

[(j) - (t)Bf- To ensure that V{4),T) > the 



Q 



A/3a vi V2 ' 



whereas the slow-rolls parameters are given by 
2 „ 2m 



e = r] 



^b) 



2 ' 



/3 



T^[l-2logi6/T)] 



^b)' 



11 



Vl 



(55) 



c = n . (56) 



Panel in Fig. [5] 
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N 


r 


top left (to = n = 0) 


60 


0.177 


75 


0.158 


top right (m = 2, n — 0) 


60 


0.225 


75 


0.193 



TABLE III: Results from the second example with A = 1 (very similar results, not shown, follow for A = 1/2). 
Rows from top to bottom refers to panels of Fig. [2] from left to right. 

To enforce the smallness of the three first slow-rolls parameters, the inflaton field must fulfill 
(/) > (/>ij + \/2. 

Thus, the energy density of the radiation fluid can be written as 



p^ = aicP-(^BrT'i2logi5/T)-l). 

Following the procedure of the previous example, we find a lower bound for Q, namely, 
n^'^'^^ ( (A A. ^2 (21og[VT)-l] 3 C 

and using ()57p . this relation simplifies to 

Clearly, in this case as well, viscosity increases the range of Q. 
The corresponding upper bound on ^O; is 



(57) 



(58) 



(59) 



Cn < C^^^ 

SU ^ Smaa 



3T 



(60) 



Figure [3 illustrates the dependence of the tensor-scalar ratio on the spectral index, from Eqs. 
dMD, and dMD- We have used the following values: a = 3 x 10-^°, 7 = 1.5, Co = (2/3)CmL, (^b = 2, 
S=l^k = (3/2) kmin and a = 10-i\ T = 2.5 x 10-^ Fq = 1.2 x 10-^ tq = 3.73 x 10-^ 0o = 0.3 
for m = 0,n = 0; and a = lO-^T = 1.75 x 10-^ Fq = 3.58 x 10-^ tq = 5.63 x 10-^ (/>o = 0.6 for 



771 = 2, n = 0, in Planck units 



M- 



Figure [2] hints that none of the two pair of indices is clearly preferred over the other. On the 
other hand, it shows that there is practically no difference between the choices A = 1 and A = 1/2. 



Table IIIII indicates the value of the ratio r for A = 1 and different choices of the pair of indices 
m and n when the number of e-folds is 60 and when it is 75. Very similar values (not shown) follow 
for A = 0.5. All of them can be checked with the help of Eqs. ()36p and (j39p . 

Again, a the results of both Tables show, in this example, the impact of viscosity is rather low. 
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FIG. 2: Top row of panels: Plot of the tensor-scalar ratio r as a function of the spectral index Ug, for two 
values of the A parameter in the case of example 2 (i.e., potential ((54)) ). Bottom row: Same as the top row 
but assuming no viscosity (Co = 0). In each panel the 68% and 95% confidence levels set by seven-year 
WMAP experiment are shown. The latter places severe limits on the tensor-scalar ratio [21]. 



Panel in Fig. [2] 
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bottom left (m = n = 0) 


60 


0.175 


75 


0.159 


bottom right (m = 2, n = 0) 


60 


0.224 


75 


0.193 



TABLE IV: Results from the second example in the absence of viscosity, i.e., (q — 0. 



VI. CONCLUSIONS 



We studied various aspects of warm inflationary universe models when viscosity is taken into 
account. The latter (whether physical or effective) is a very general feature in multiparticle and 
entropy producing systems and, in the context of warm inflation, it is of special significance when 
the rate of particle production and/or interaction is high. In consequence we have focused on the 
strong regime Q ^ 1 (in the weak regime viscosity is, at best, negligible [19]). 

Equation (j22p gives the necessary and sufficient condition for the existence of stable slow-roll 
solutions. Under this condition, we got the same stability range obtained in the no-viscous case, 
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so long as o" = —8/3. In this sense, the range of the slow-roll parameter c decreases when —8/3 < 
(T < 0, and increases when a < —8/3. 

We calculated a general expression for the spectral index, Ug (Eq. ([HU]) ). that depends explicitly 
on viscosity through the four pi coefficients. The latter do not depend on the slow-roll parameters 
(e,/3,r/, and6), as shown by equations ([3T|) - ([M|) . 

In order to further ensure the stability of the warm viscous inflation, the slow-roll parameters 
must satisfy the following conditions 

e<Cl + (3, |/3|<1 + Q, \r]\<^l + Q, 

as well as the condition on the slow-roll parameter that describes the temperature dependence of 
the potential, namely, 

To get explicit expressions we have considered two examples for the scalar potential. In the first 
case, V{(j),T) = vi{(j)) + V2{T) (Eq. (j42p ). it trivially follows that b vanishes; on the other hand, 
the stability condition requires that the strength parameter fulfils 

Q > 9*-Ps [I - -^] . 

\ smax J 

Hence, the viscous pressure increases the stability range of Q as Co < Cmax, where Qnax is given by 



Eq. dMl). 

The second example is provided by the potential V[4>^T) = avi{(j))v2{T), (Eq. ([M|) ). In this 
case, the stability condition for Q is 



^ > 42^4 I ^ (2) ) ' 

\ '^max / 



with (max given by ([60]). Again, the range of stability of Q augments in the presence of viscosity. 

To bring in some explicit results we have taken the constraint Ug — r plane to first-order in the 
slow roll approximation. For the first potential (Eq. (j42p ) we obtained that, when A = 0.5 and 
A = 1, the model is consistent with the WMAP seven year data for the pair of indices {m = 2, 
n = 0), see Fig. [TJ Likewise, for the second potential (Eq. (j54p ) all two choices of the (m, n) pair 
show compatibility with the said data -see Fig(TJ The results obtained indicate that the effect of 
viscosity on the value of r is marginal. However, it cannot be discarded that future experiments 
uncover it. 

We should note that there are some properties in this model that deserve further study. In 
particular, we did not address the case n ^ 0. In fact, when F = F(0, T) a more detailed and 
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laborious calculation for the density perturbation would be necessary in order to check the validity 
of expression (j28p . We intend to return to this point in the near future by following an approach 
analogous to that of Ref . [201] . 
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